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Abstract. A 180 /im map and strip maps at 120 and 180 fj,m were obtained for the edge-on starburst galaxy NGC 
253 with ISOPHOT, the photometer on board the Infrared Space Observatory. We compare these observations 
with those obtained by IRAS at 60 fim and 100 /im and derive the far-infrared spectral energy distribution at 
different locations in the galaxy. There is evidence for the presence of cold dust (T < 20 K) both in the nucleus 
and in the disk. Extended emission dominated by cold dust is detected up to ~ 15' (~ 10 kpc) along the major and 
minor axis; its spatial distribution is similar to that seen in the IRAS and ROSAT PSPC images. The emission 
along the minor axis is probably related to large-scale outflows of gas (superwinds) which originate in the nuclear 
starburst and maybe to star formation in the halo. The radial dependence of the dust temperature along the 
major axis is found using a radiative transfer code: we show that the dust scale length in the disk is ~ 40% larger 
than that of stars. 
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1. Introduction 



There is growing evidence (Trewhella et al. 200Cl| ) for the 
existence of two components in the dust distribution in 
galaxies, a spatially extended distribution of cold dust 
(T < 20K) and a warm dust (T > 20K) component con- 
centrated in the disk, probably associated with star for- 
mation. The most direct evidence for the existence of 
the cold dust component has been provided by SCUBA 
and by ISOPHOT, the photometer on board the Infrared 
Space Observatory (ISO, Kessler et al. 1996), in particu- 
lar thanks to the C200 detector. The cold dust may lie 
both along the disk, in which case the dust scale length 
may be higher than the stellar scale length by factors of 



40% to 80% (Ahon et al. |1998| b) and vertically. If a mas- 
sive nuclear starburst is present, supernovae forming there 
may drive the gas from the inner disk of the host galax- 
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iesto their halos, as first propose d by C hevalier & Clegg 
( 1985 ) and later Heckman et al. ( 1990 ). Dust entrained 
in these outflows should enhance the far-infrared emission 
vertically to the disk; this effect should be most easily 
observed in edge-on galaxies. 

NGC 253 is a nearby, almost edge-on barred spiral 
galaxy (De Vaucouleurs et al. 



1991); we adopt a distance 



of 2.5 Mpc (Houghton et al. |1997|) , which gives a projected 
linear scale of 12pc/arcsec. NGC 253 is one of the nearest 
galaxies showing evidence for a compact nu clear starburst 
(size ~ 100 pc, Dudley & Wynn- Williams 1999 ) which is 
partially responsible for its high far-infrared luminosity, 
Lfir ~ 2 X IO^'^Lq. NGC 253 has been extensively ob- 
served in all spectral r anges (see e.g. Engelbracht et al. 
1998 and Forbes et al. 200C ). Still uncertain is the pres- 



ence of a highly obscured AGN which would contribute 
some of the radio and X-ray emission, but none of the 
optical and infrared emission (Forbes et al. 2000). 



McCarthy et al. ( 1987 ) found along the minor axis an 
extended emission-line region associated with diffuse x- 
ray emission, over ~ 15" in diameter: they identified this 
emission-line gas with cool dense clouds embedded in an 
outflowing wind of hot gas. Fabbiano (1988) found more 
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extended, low surface brightness X-ray emission extending 
along the minor axis out to 10' (NW) and 5' (SE). More 
recently, Dahlem et al. ( |1998D analyzed ROSAT PSPC 
and HRI and ASCA data and found a diffuse, hot soft 
X-ray halo with dimensions of 16 x 10 kpc and attributed 
it to the superwind. Pietsch et al. (2000) also investi- 
gated the extended X-ray emission from ROSAT HRI and 
PSPC data: they separated the total halo emission in a 
corona component which originates from the halo imme- 
diately above the disk (scale height ~ 1 kpc) and in an 
outer halo which extends to projected distances from the 
disk of 9 kpc and shows a horn-like structure. However, 
IRAS images of NG C 253 did not provi de stra ightforward 
results. Rice ( 1993 ) and Alton ct al. (1998i) presented 
high-resolution IRAS maps produced with the Maximum 
Correlation Method (MCM) technique; they concluded 
that spurious effects like detector hysteresis at 12 and 
25 fim and emission reflected from the IRAS telescope spi- 
der at 60 and 100 /im mask low-level emission outside the 
disk. Alton et al. (1999) presented submillimeter images 
of the nuclear regions of NGC 253 obtained with SCUBA 
and found evidence for a dust outflow along the minor axis 
at 450 /im on a size < 45". 

In this paper we report the results of observations of 
NGC 253 obtained with ISOPHOT; available data and 
their processing are discussed in sec. |[ The investigation 
aims at: 

1. Evaluating the far- infrared emission properties in dif- 
ferent regions of the galaxy and look for evidence of 
the outflow in the FIR. To this purpose we combine the 
ISOPHOT 180 Aim map with the IRAS maps (sec. |): 
an estimate of dust temperatures and masses is ob- 
tained by fitting the so-obtained spectral energy dis- 
tributions with a simple model (sum of two modified 
blackbodies). 

2. Analyzing the radial dependence of the dust tempera- 
ture using a radiative transfer code (sec. FIR emis- 
sion properties in the disk are compared to those de- 
rived from optical data, which are mainly related to 
the stellar component. 



2. Observations and Data Reduction 

NGC 253 was observed with ISOPHOT (Lemke et al. 
1996| ) in the staring raster mode with the multi-filter 



Astronomical Observation Template PHT22 (Klaas et al. 



1994) and the C200 detector. The staring raster mode is 



a sequence of staring observations on a two-dimensional 
(M X N) regular grid which consists of a sequence of space- 
craft pointings. As given in Table ^ a map was obtained 
at 180 /zm covering a 53' x 26' field. The separation of 
the raster points was equal to the C200 array size (180"), 
which resulted in gaps between them. In addition, four 
strip (N=l) maps were performed in the 120 and 180 /im 
filters. They were oriented roughly along the minor axis 
of NGC 253 and centered on the nucleus and offset posi- 
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Fig. 1. The areas covered by ISOPHOT maps have been 
overlaid to J band image of NGC 253 from the Digitized 
Sky Survey, obtained using the SkyView survey analysis 
system. 



tions of 10, 15 and 20'. The positions of the four maps are 
shown in Fig. |l|. 

The data reduction was performed using the 
ISOPHOT Interactive Analysis tool (PIA, version 7.2) to- 
gether with the calibration data set V 4.0 (Laureijs et al. 
1998): corrections were made for non-linearity effects of 
the electronics, disturbances by cosmic rays (deglitching) 
and signal dependence on the reset interval time. The flux 
calibration is based on measurements with the thermal 
fine calibration sources (FCS) on board: two FCS mea- 
surements were taken, one before and one after each ob- 
servation in order to correct for changes in the detector 
responsivity According to Klaas et al. ( |2000[ ), ISOPHOT 
absolute calibration accuracies for extended sources ob- 
served in the staring raster mode are < 25%. Saturation 
occurred in the nucleus both at 120 and 180 /zm; this could 
be at least partially recovered only in the 180 /zm strip map 
including the destructive read-outs, which are usually dis- 
carded. 

Our data were further processed using a series of pro- 
cedures developed by Manfred Stickel at MPIA. A first set 
of procedures allowed to improve the flat field correction 
of the 180 /im map. A second set of procedures allowed to 
interpolate between the detector gaps changing the pixel 
size to 50" and to rotate the map; these procedures are 



based on the 'drizzle' algorithm (Fruchter & Hook 1997) 
implemented in IRAF. ISOPHOT maps of a point-like 
source (NGC 7027) at 120 and 200 /im were processed fol- 
lowing the same procedure, in order to measure the size 
(FWHM) of the beam profile: in both cases we obtain 
FWHM = 122". 
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Table 1. Log of the ISOPHOT observations. The position angle (PA) is from North to East. Offsets from the nucleus 
are given for the strip maps. 
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Fig. 2. Surface brightness profiles normalized to the peak brightness in the strip maps at 120 /zm (open circles) 
and 180 /xm (filled squares) at offsets of 0', 10' (left) and 15' {right). The lower limit has been chosen to emphasize 
the extended part of the emission; error bars are computed from the statistical uncertainties on measured surface 
brightnesses. The dashed line is the beam profile derived from ISOPHOT maps of a point-like source at 120 and 
200 ^m. 



2.1. IRAS maps 



3. Mapping of the FIR emission 



IRAS images of NGC 253 processed with the High 
Resolution (HIRES) technique were obtained from the 
Infrared Processing and Analysis Center (IPAC). A de- 
tailed description of this technique was given by Rice 
( 1993 ). We used the maps with twenty iterations, giv- 
ing resolutions of 37x23" (12 ^m), 35x23" (25 /im), 62x41" 
{eOfim) and 98x80" (100 ^im). 

The IRAS maps were convolved with a Gaussian func- 
tion to the FWHM of the ISOPHOT beam profile at 
200 Atm (122") and registered to the ISOPHOT 180 ^m 
map using the coordinate information contained in the 
headers of the image files: the registration also included 
matching the pixel size in the maps. The reliability of the 
registration was checked by the coincidence of the surface 
brightness peak; the IRAS and ISOPHOT maps are com- 
pared in Fig. ||. 



The surface brightness profiles in the strip maps show (see 
Fig. ^ that the emission is dominated by the strong nu- 
clear source. However, extended emission is also present 
along the minor axis, in particular in the strip centered 
on the nucleus (0' offset) where it may be detected within 
uncertainties to a distance ^ 8' from the nucleus and with 
decreasing signal to ~ 12'. We define this emission as the 
'halo' component, in addition to that detected along the 
major axis ('disk') and in the center ('nucleus'). Along 
the major axis, emission is detected up to 15', where the 
data at 120 /im are significantly noisier than at 180 /zm; no 
emission above the background is detected in the 20' strip 
map. In order to measure the fluxes in these three compo- 
nents, we defined circular apertures at different positions 
in the galaxy (see Fig. || and Table ||) and computed the 
integrated fluxes from the background-subtracted IRAS 
and ISOPHOT maps. 

The flux measured in the aperture centered on the nu- 
cleus is the sum of disk and nuclear emission. Due to the 
poor spatial resolution it was not possible to separate the 



4 



M. Radovich et al.: Far- infrared mapping of NGC 253 




30*^ 0''48'"00'^ 30^ o''47™00'^ 30^ 30^ o''48™00'^ 3 0*'47™00'^ 30^ 

a (2000.0) a (2000.0) 



Fig. 3. Contour maps of the IRAS and ISOPHOT images and circular apertures where fluxes have been measured. 
The IRAS maps were transformed to the same position and resolution as the ISOPHOT map. The contour levels 
(MJy/sr) are: 12 ^m - 0.8, 1.5, 2.5, 5, 10, 30, 50, 80, 300, 500, 1000, 1500; 25 ^m - 1.2, 1.5, 2.5, 5, 10, 30, 50, 80, 300, 
500, 1000, 1500; 60 /im - 1.0, 1.5, 2.5, 5, 10, 30, 50, 80, 300, 500, 1000, 1500; 100 ^im - 1.5, 2.5, 5, 10, 30, 50, 80, 300, 
500, 1000, 1500; 180 /im - 1.5, 2.5, 5, 10, 30, 50, 80, 200, 300, 370. 



two components in our data. We made an attempt to esti- 
mate the emission that would be produced by a 'point-like' 
source (FWHM < 122"); we shall define this component 
as the 'nucleus', with the understanding that actually it is 
the sum of nuclear and disk emission inside the region that 
would be seen as 'point-like' by the C200 detector. To this 
aim, we follo wed a procedure similar to that described by 
Ahon et al. (|l998| a). We normalized the 2D beam profiles 
to the observed peak surface brightnesses and integrated 
them in an aperture of diameter ~ 8'. In the case of the 
IRAS images we used the beam maps provided by IPAC. 
In the case of ISOPHOT we used a 2D gaussian fit of the 
point-like source where the beam profile was measured: 
this was normalized to the peak surface brightness found 
in the 180 /im strip map since the nucleus in the large 
180 ^m map is saturated. The ratio of the IRAS total to 
nuclear (rescaled beam maps) fluxes gives results similar 



to those found by Alton et al. (1998a) for infrared-bright 
galaxies, namely -FeoV^eo"'' ^ 1-4 and Fl°yFl]^^ - 1.8; 
from the ISOPHOT map we get - 3. It follows 

that while ^ 70% of the 60 fim emission comes from the 
nucleus, emission at longer wavelengths is dominated by 
the outer regions. 

In the halo component, fluxes from the IRAS maps 
may be affected in an unknown way by instrumen- 
tal spurious effects, as seen above. As explained in 
www.ipac.caltech.edu / ipac / iras / iras_data_features.html, 
very bright (> 500 Jy) pointlike sources may produce 
a six-pointed star shape due to the reflection from the 
secondary mirror struts; there is no known method to 
remove this effect. Such a pattern is seen in the IRAS 
maps of NGC 253 at 60 ^m and 100 /im (Fig. |^); however, 
we note that the same pattern in the 'putative' diffuse 
emission may be recognized not only in the ISOPHOT 
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Fig. 4. Spectral energy distributions measured in the regions displayed in Fig. H; fluxes have been color corrected. 
IRAS fluxes are displayed by squares, ISOPHOT fluxes by circles; crosses show ground-based phot ometry of the 
nucleus (5.5" aperture) taken from NED; the triangle is the 1300 /Ltm flux given by Kriigel et al. (1990). 



180/xm map, but also in the ROSAT images (see e.g.: 
Pietsch et al. ^GOOt panels B (0.1-2.4 k eV) a nd S (0.1-0.4 
keV) in their Fig.l; Forbes et al. ^OOOj Fig.2). We 
therefore conclude that the diffuse emission detected 
in the IRAS maps is probably real at A > 60 /im and 
maybe at 25 /im, while no significant contribution is 
detected at 12 /im; since we cannot exclude the existence 
of instrumental effects, fluxes from low surface brightness 
regions in the IRAS maps should be considered as upper 
limits. 



3.1. Spectral Energy Distributions 

Physical quantities (dust temperatures and masses) in the 
different regions were first derived fitting the observed 
SED with modified blackbodies, i.e. blackbodies with 
a wavelength dependent emissivity (see e.g. Hildebrand 



1983). Far-infrared fluxes in the obs erved range (12- 
180 /im) are probably due (Desert et al. 199Q|) to a mix of 



polyaromatic hydrocarbon molecules (PAHs), very small 
grains and big grains. Emission is dominated by PAHs at 
12 /im, by co ld sili cate and grap hite g rains at A > 140 /im 
(Dwek et al. |l997l Sturm et al. |2000| , Genzel & Cesarsky 
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Table 2. Distances from nucleus (r) and diameters (d) of the apertures within which fluxes were measured; dust 
temperatures, infrared luminosities (1-1000 /im) and dust masses derived from the fits with and without the 25 /um 
point are displayed. 



Total 


Nucleus 


Disk 


Halo 


a 


b 


c 


d 


e 


f 


rn 





8 


8 


9 


8 


r [kpc] 





6 


6 


7 


6 


d ['] 33 


8 


5 


5 


8 


5 


d [kpc] 24 


6 


4 


4 


6 


4 


Observed fluxes, not color corrected 










12 /im 63.2 


35.0 


2.0 


1.7 


- 


- 


25 /im 181 


178 


2.4 


1.5 


3.3 


< 0.1 


60 /im 1186 


853 


19.1 


12.2 


6.1 


3.6 


100 fim 1770 


949 


56.6 


44.0 


25.5 


8.3 


180 fim 1693 


607 


87.3 


77.4 


28.3 


15.7 


Fit from 25/1171 to 180 fi 


m 










Warm component 












T [K] 60 


60 


60 


60 


80 


40 


Lir [Lq] 1.07E-I-10 


8.18E-H09 


1.44E-h08 


9.08E+07 


8.32E-f07 


3.33E-f07 


Md [M0] 7.42E-h04 


3.59E+04 


1.04E-h03 


6.47E+02 


6.75E-f01 


1.61E+03 


Cold component 












T [K] 20 


20 


20 


20 


20 


15 


Li, [Lq] 8.51E-f09 


4.24E-I-09 


3.76E+08 


3.14E+08 


1.56E-h08 


5.80E-)-07 


Md [Mq] 3.51E+07 


7.23E+06 


2.69E+06 


2.79E+06 


5.42E-f05 


9.84E+05 


Fit from 60 jim to 180 


m 










Warm component 












T [K] 47 


39 


28 


27 


25 


29 


Lir [Lq] 7.90E-h09 


5.76E-H09 


1.91E-h08 


1.27E-h08 


1.26E-H08 


4.60E-H07 


Md [Mq] 2.36E-h05 


5.47E-f-05 


1.34E+05 


1.05E+05 


1.86E+05 


2.32E-f-04 


Cold component 












T [K] 20 


20 


17 


16 


17 


13 


Lir [Lq] 7.72E+09 


2.22E-H09 


2.51E-h08 


2.29E-h08 


4.40E-f07 


4.22E-f07 


Md [Mq] 3.83E-I-07 


9.81E-f06 


3.88E-h06 


3.85E+06 


6.74E+05 


3.30E+06 



200C ) ; at intermediate wavelengths each of the three com- 
ponents may contribute to the observed fluxes. Our pur- 
pose is to determine temperatures and masses of the cold 
dust component, but due to the limited number of points 
it is not possible to obtain a univoque fit (see Lisenfeld 
et al. 20GC for a discussion about the non uniqueness of 
single-temperature fits for the far-infrared SEDs). In order 
to estimate the impact of this uncertainty we first made a 
fit including all the points with the exception of the one 
at 12 /im, which is dominated by PAHs. The fit was then 
repeated dropping the 25 /xm point, leading to lower dust 
temperatures and higher dust masses; values from this fit 
are indicated in italics. 

We used for the fit the sum of two modified blackbod- 
ies and an emissivity oc A^^. Color correction was applied 
according to the temperatures found and the whole proce- 
dure was repeated until convergence. The observed fluxes 
and fitted SEDs are shown in Table ^ and Fig. |^. For each 
component, we computed dust masses according to Klaas 
& Elsasser (|1993D , assuming dust grain properties as in 
Hildebrand ( |1983D for the emissivity adopted here : 



where Lir is the luminosity between 1 and 1000 /im 
computed from the extrapolation of the modified black- 
bodies. 

Total fluxes from the whole galaxy (a) were computed 
from the integration in a circular aperture with 33' diam- 
eter. Fluxes from 25 to 180 /tm are well fitted by a two- 
component blackbody with T ~ 60(.^?) K and T - 20 
K. Using the dust mass derived from the SED integrated 
over the whole galaxy (Md ^ 3.5(5.<S) x 10^ Mq) and the 
total mass (Af(HI + H2) = 2.4 x 10^ Mq) estimated by 



7.9 X 10-^(Tk/40)-^ im/io [Mq], 



(1) 



Houghton et al. (1997) we obtain a dust to gas mass ratio 
^ 1/70, somewhat higher than the Galactic value (1/160). 

1. Nucleus - The nuclear SED (6) is fitted by two com- 
ponents, a warm one (T ~ 60(5P) K) and a cold one 
(T ~ 20 K). The ratio of the extrapolated 1-1000 /im 
luminosities in the two components is L^/Lc ~ 2, i.e. 
the emission is dominated by the warm component. 
We estimate a dust mass ~ 7(JO)xlO^ M©. In order 
to check the reliability of the separation of the 'nuclear' 
emission from the extended emission, we extrapolated 
the nuclear SED to 1300 /i m and compared it with the 
flux given by Kriigel et al. ( 1990 ), 1.7 Jy in an area of 
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3' X 2': the two values agree within 30 (<§)%, which is 
acceptable considering the uncertainties involved. 

2. Disk -The SEDs in the NE (c) and SW (d) regions of 
the disk show, as in the nucleus, the presence of a warm 
(T ~ 60(28) K) and a cold (T ~ 20(i7) K) component. 
However, emission is dominated by the cold compo- 
nent, L^/Lc ^ 0.3(0.5). The dust mass estimated in 
each of these regions is ~ 3(-^)xl0® M0. 

3. Halo - Due to the uncertainty in the IRAS fluxes, a 
SED decomposition is very difficult; the results given 
here are therefore indicative. In the southern halo (/) 
we could only assume an upper limit for the flux 
at 25 /im since no significant emission was detected. 
Dust temperatures of the cold component are simi- 
lar to those found in the disk, T ^ 20(-Z7) K, with a 
somewhat lower temperature in the southern halo (/) 
where T ^ 15(^5) K; we derive dust masses in the or- 
der of 0.5-1(0. 7-5) X ID*' Mq. In order to explain the 
non detection by IRAS of emission by dust involved 
in outflows, Alton et al. ( |l998| a) argued that the tem- 
perature of the dust is T ~ 15 — 20 K and its mass 
< (0.5 - 3.5) X 10^ Mq: this hmit is not far from the 
values we find for NGC 253, in particular if we consider 
that if T 20K an uncertainty in the temperature of 
AT = 5K may change the dust mass by a factor of 
5. For comparison, Phillips (1993) gave a dust mass of 
2 X 10 ^ Mq for the outflow in NGC 1808; Ahon et al. 
( 1999j ) estimated for the outflow in M82 dust masses 

~ 10^ M„. 



of 10'' 



According to Kennicutt (1998), the star formation rate 
for continuous bursts of age 10-100 Myr in starbursts, 
where contribution of dust heating from old stars may 
be neglected, is: 



SFR = 1.7 X 10-i%r/Lq [Mq yr-i]. 



(2) 



The SFR obtained in the nucleus and in the whole galaxy 
is ~ 2.1(i.5) and 3.3(^.^) Mq yr-\ respectively. This low 
value of the SFR i s in a greement with the conclusion by 
Engelbracht et al. (|l998| ) that NGC 253 is in a late phase 
of the starburst, having passed a rapid decrease of the star 
formation rate. 



4. Modeling the FIR emission 

The far-infrared emission of NGC 253 was modeled with a 
Monte Carlo simulation system (Kahanpaa 2001, in prepa- 
ration) using both the 180 /zm map for FIR major and mi- 
nor axis profiles and the strip maps at 0, 10 and 15' for 
minor axis profiles and photometry at 120 and 180 /im. 



4.1. The galaxy model 

The galaxy model is based on a cylindrically symmetrical 
disk with exponential distributions along the radius and 
height with an additional, compact nuclear energy source. 
The distribution of stars within the disk was assumed to 



be independent of wavelength; no intrinsic color gradi- 
ents in the stellar population were taken into account. 
While this ignores color gradients due to the bulge, the 
impact on the resulting FIR colors is minor as NGC 253 
is a disk-dominated spiral (Baggett et al. 1998). The core 
was treated as an independent source with its own unob- 
scured spectral energy distribution. Since recent studies 
by Davies et al. (1999) suggest that the scale lengths of 
star and dust populations are not equal, the scale lengths 
were treated as independent parameters in the model. 

4.2. Radiative transfer code 

Unlike most radiative transfer codes our model does not 
ignore the cloud structure of the interstellar medium 
(ISM). Instead, all dust is in clouds; the total number of 
clouds is determined by the total dust mass and the cloud 
positions are randomly chosen from a double-exponential 
distribution: 

/J'Me-"/''«du 

-Rl = rr„„^ -, ; (3) 



i?2 = 



3-"/''odw' 



(4) 



where Ri and R2 are uniform random numbers and tq and 
zq the dust scale length and scale height respectively. 
The radiative transfer code for the optical domain is 



based on models by de Jong (1996) and Bianchi et al. 
( 19961 ). The assumed unobscured spectral energy distri- 
bution of the energy sources (stars) was divided into a 
number of energy bins - each formally representing one 
broad-band filter - and the radiative transfer problem 
was numerically solved for each bin by following numerous 
simulated photons through the dust system. For each bin 
2 • 10^ to 1 • 10^ random paths through the dusty layer were 
calculated; the total number of independent paths used in 
a single model was typically 40 • 10^. The resulting surface 
brightness distribution for each wavelength was recorded. 
Additionally the cumulative amount of absorbed energy 
at each position within the model galaxy was registered 
for further processing. When all bins had been processed, 
the equilibrium temperature of individual dust grains as a 
function of grain size and position within the galaxy (dis- 
tance from the nucleus and height above the plane of the 
galaxy) was solved. In thermal equilibrium the emitted 
power from a single grain is 

/>OC 

Pem= / 4^'a2Q,bs(A,a)B(A,T)dA, (5) 
Jo 

where a in the grain radius, Qabs is the absorption coeffi- 
cient of the grain and B is the blackbody (Planck) func- 
tion. The shape of Qabs in the FIR was approximated with 

Qabs(a,A) = 1.8 X lO-^aA-'^. (6) 

We have used /3 = 2, a typical value for grain tem- 
peratures less than 50 K, and dust grain size distribu- 
tions and optical properties from Desert et al. (1990). 
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radius (kpc) 
4 6 




Table 3. Model input parameters and results for NGC 
253 



Fig. 5. Grain temperature in the galactic plane of NGC 
253. The solid line represents grains with a diameter of 
150 nm, the dotted line the smallest classical grains in the 
model (d = 15 nm). 



The dust model chosen is of secondary importance; we 
tried using the classical MRN dust population models with 
n{a) oc a~^-^ (Mathis et al. 1977) and got the same results 



for classical dust emission. Infrared emission from very 
small grains and PAHs with stochastic heating and cooling 
processes were neglected as they are minor contributors in 
the 100-200 /xm wavelength range studied here; they were, 
however, still included as absorbers with size distributions 



and optical properties as given by Desert et al. ( 1990 ) 



4.3. Model input parameters 

Table || summarizes the parameter values used in the 
model. The observed total stellar fluxes as a function of 
wavelength were obtained from the NASA-IPAC extra- 
galactic database (NED). The radial scale length of the 
stellar component was obtained by comparing the result- 
ing brightness profiles for the optically thin outer parts of 
the axial profiles with published measurements of the ma- 
jor and minor axis profiles in the B and V bands (Pence 



198C, Baggett et al. 1998). As a byproduct the inclination 
of NCG 253 was found to be 78°± 1°. Limiting values 
for the radial dust scale length were found by comparing 
the model results with the observed major axis profile at 
180 /im. Since the model profile depends also on the tem- 
perature of the dust grains, this fit does not give as strin- 
gent limits as in the optical case, but if the scale length 
was varied by more than 15%, the fit was unacceptably 
poor. The major parameters of the model were adjusted 
until a reasonable fit to the optical and FIR data was 
found. Unobscured SEDs for the disk and nucleus com- 
ponents were appr oxima ted with model SEDs for stellar 
populations (Silva 1998). For the core we used a 200-My 
old population without modifications; in the disk a com- 
bination of young (0.2 Gy) and old (5 Gy) populations 



Input pnvanijCtcT'S 




Distance 


2.5 Mpc 


Stars, luminosity 


3.2 X 10^° Lq 


Stars, core to disk luminosity ratio 


0.5 


Dust, cloud size 


40 pc 


Dust, grain density 


3 gcm"^ 


Dust, maximum albedo 


0.7 


Dust, size distribution 


As in Desert et al. 1990| 


Results 




Inclination 


78° 


Stars, scale length 


1.8 kpc 


Stars, scale height 


0.2 kpc 


Dust, scale length 


2.8 kpc 


Dust, scale height 


0.15 kpc 


Dust, total mass 


2.8 X 1O''M0 


TV (center) 


4.5 ± 0.5 


FIR luminosity 


1.7 X 10^" Lq 


Large grain Lumin. 


1.2 X 10^" Lq 


VSG & PAH Lumin. 


0.45 X 10^" Lq 


T(r=0.5 kpc) 


25 ± 1 K 


dT/dr in disk 


1.2 Kkpc"^ (1-10 kpc) 



was required. The resulting unobscured SED is shown in 
Fig. I 

4.4. Predicted dust mass, distribution and temperature 
profiles 

According to the computed model the face-on optical 
depth of the dust layer at the center of the galaxy is 
4.5 at 550 nm. After the inclination of NGC 253 (78°) 
is taken into account this makes the core attenuated by 
roughly 11 magnitudes in the V band. Values based on 
earlier observations are somewhat higher: Engelbracht et 



al. ( |1998D estimated Av(cen) = 19.07 ± 3.3 from the H-K 
excess, but as the small starburst core certainly includes a 
sizeable amount of hot local and circumstellar dust, this is 
to be expected. Half (54%) of the total model luminosity 
(3.2 X 10^° Lq) is converted to FIR emission by the dust 
grains. The total mass of the cold dust component as in- 
ferred from the model is M ~ 0.93 x lO^pd Mq, where pd 
is the average density of the dust grains in units of g cm~^ . 
Using the value (pd — 3 gcm~^) given by Hildebrand 



( |1983D sets the dust mass of the model to 2.8 x IO^Mq, a 
value not too far from the estimate found by using Eq. |[ 
A disk-to-nucleus luminosity ratio of 2.0 gave the best 
match with the observed SED and also results in a core 
luminosity Lc 



1.1 X IO^^Lq; the same value was found 



by Engelbracht et al. ( 199^ ) from a combination of mid- 
and far-infrared, sub-millimeter and 1.3 mm data. 

Fig. ^ shows the radial temperature distribution of 
dust in the plane of NGC 253. The smallest grains tend to 
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Fig. 6. Spectral energy distribution of NGC 253 in 
the optical-near infrared (ground-based photometry from 
NED) and far-infrared (IRAS 100 ^m and ISOPHOT 
180 ^m fluxes). Crosses - observed SED; squares - unob- 
scured SED; dashed line - model fit for the emission from 
stars and cold dust, as described in Sect. |[ IRAS fluxes 
from 12 to 60 /im were ignored in the fit since the model 
does not include the IR emission from PAHs or very small 
grains, which are responsible for most of the emission at 
these wavelengths. 



be hotter and thus more efficient FIR emitters than their 
bigger companions because the FIR absorption coefficient 
of the grains at a given wavel ength is proportional to grain 
diameter (see Draine & Lee 1984 ). This effect is partially 
counterbalanced by the spectral shape of the attenuated 
interstellar radiation field (ISRF) within the galaxy: the 
smallest grains are poor absorbers in the near infrared, 
where the ISRF is most intense. 

The central dust temperature given by the model of 
NGC 253 is higher than that given by a simple fit to the 
FIR spectrum (Fig. ^). The central spectrum is, however, 
integrated over quite a large area and in the end the fit is 
sufficiently good. The hot central core is very small and 
its contribution to the total SED is approximately 40% of 
the total luminosity. The total SED (Fig. ||) as well as the 
strip map flux levels at 120 and 180 /im at 0' and 10' are 
perfectly reproduced. 

The derived scale lengths of the stellar and dust com- 
ponents are 1.8 and 2.8 kiloparsecs respectively. As the 
uncertainty of both values is of the order of 15%, the dif- 
ference is significant. This agrees with the recent observa- 
tions of dust dist ribution in ot her n earby spiral galaxies 
(Alton et al. |l998| b, Alton et al. |200l| ). A prominent effect 
of the difference in the scale lengths is the variation of 
the dust to stellar mass ratio with radius: at 10 kpc the 
relative dust density is four times as large as at the center. 

The halo dust component cannot be reproduced by a 
single 3-D exponential dust disk: the modeled minor axis 
profiles (convolved to 120" FWHM) are narrower than the 
observed proffies of NCG 253 at 120 and 180 /xm. 



This can be due to several different effects or a combi- 
nation of them: 

1. There is an extra source heating the halo, which is not 
represented by the stellar populations in the model. 

2. A halo population of young stellar objects (YSO) or 
other compact FIR emitters is present in NGC 253. 

3. The dust distribution in the halo has a distribution 
which significantly deviates from a simple exponential 
model. 

The two first scenarios are in agreement with the exis- 
tence of star formation in the halo of NGC 253, which has 



been recently proposed by Comeron et al. (2001) on the 
basis of UBV photometry. In this case we should be able 
to see the optical emission of the heating sources as the 
halo is certainly optically thin. The very cold spectrum of 
the outer parts contradicts the second case as YSO sources 
usually have SED with strong near and mid-infrared com- 
ponents. In the third case, which we find most probable, 
the excess dust in the halo is heated by stars in the disk 
of NCG 253. If a dust halo is the most important contrib- 
utor it should also be observable in the optical domain as 
a faint, blue halo of scattered light. 

5. Conclusions 

ISOPHOT observations were used together with IRAS 
data to map the far-infrared emission in the starburst 
galaxy NGC 253. The resulting images were then analyzed 
with modified blackbody fits as well as radiative transfer 
modeling. The total FIR luminosity consists of two parts: 
a warm component (T ~ 47-60 K) , dominated by the cen- 
tral regions, and a cold component (T ~ 20 K) produced in 
the disk and the halo. Dust temperatures and masses de- 
rived from single-temperature fits of the SED integrated 
over the whole galaxy are in agreement with those ob- 
tained by solving the radiative transfer equations. The re- 
sults obtained in the different regions may be summarized 
as follows: 

Nucleus: the nucleus produces half of the total FIR 
luminosity of NGC 253. Our fit for the core correctly pre- 
dicts the observed 1.3 mm continuum emission given by 
Kriigel et al. (1990). The low inferred star formation rates 
suggest that the nuclear starburst in NGC 253 is in a late 
phase, having passed a rapid decrease of the star forma- 
tion rate. 

Disk: the resulting dust temperatures are typical for 
the general diffuse ISM in spiral galaxies. The radiative 
transfer model indicates that the dust scale length is ^ 
40% larger than that of stars; the same trend was found 
by Ahon et al. ( |1998| b). The dust is warmer {T > 20K) 
in the central regions of the galaxy, the outer regions are 
dominated by colder dust (T < 15K). 

Halo: we find evidence for extended emission to pro- 
jected distances of ~ 10 kpc from the nucleus along the 
minor axis both in the IRAS and ISOPHOT maps; the 
emission seen in the IRAS maps could be due to instru- 
mental effects, but since the same emission pattern is also 
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seen in the ISOPHOT map as well as in ROSAT images 
(e.g. Pietsch ct al. |2000| , Forbes et al. |200C| ) we conclude 
that this is not the case. The SED fitting of this extended 
component is very uncertain, since IRAS fluxes may still 
he affected by instrumental effects. Keeping in mind this 
uncertainty, the emission is probably dominated by cold 
dust (T ~ 15 — 20 K); the derived dust masses are com- 
parable to those found for the outflows in NGC 1808 and 
M 82. The apparently constant mass in the outflows ob- 
served in nearby galaxies is probably due to limitations 
in the instrumental resolution and sensitivity. It was not 
possible to reproduce the halo component by a single 3-D 
exponential dust disk. Even if in our opinion this is most 
likely due to a non exponential distribution of dust in the 
halo, some extra heating may be give n by star formation, 
which according to Comeron et al. ( 2001 ) may occur in 
the halo. Both phenomena may be related to the interac- 
tion of the interstellar medium with the superwind driven 
by the nuclear starburst. 
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